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Abstract 

o 

5^ As a liquid is progressively supercooled an intriguing weakening of the wavenumber (q) dependence of the structural 
relaxation time r(q) in the large q limit is observed both in experiments and simulation studies. Neither the continuous 
Brownian diffusive dynamics nor the discontinuous activated events can alone explain the anomalous wavenumber 
dependence. Here we use our recently developed theory that unifies the mode coupling theory (MCT) for continuous 

^ dynamics and the random first order transition theory (RFOT) treatment of activated discontinuous motion as a 
nucleation (instanton) process, to understand the wavenumber dependence of density relaxation. The predicted 

s 

smooth change in mechanism of relaxation from diffusive to activated, in the crossover regime, is wavevector dependent 
and is eventually responsible for sub-quadratic, almost linear, q dependence of the relaxation time. 

o 

I. INTRODUCTION 

(N 
> 

oo . . . 

£■ — Diffusion in an equilibrated supercooled liquid, whatever the mechanism of transport at the microscopic level, should 

o 

follow Fick's law at long length scales and therefore, the relaxation time should exhibit a quadratic q dependence at 

(N 

small q. However, relaxation at small (or molecular) length scales (large q) is sensitive to the detailed mechanism of 

transport. In deeply supercooled liquid, there is a crossover from continuous (in space) small amplitude Brownian 
> r 

. movement of particles to rare, nearly discontinuous, larger amplitude hopping motions. This crossover appears to 

X 

occur over a range of temperature which is limited on one side by the temperature (Tl) that signifies the onset of 
barriers for motion in the inherent structures, and ultimately by the glass transition temperature on the other side 
when the system falls out of equilibrium [1 . It has been shown that this onset temperature Tl [2] coincides with 
the ideal mode coupling theory (MCT) temperature T° [3]. In Salol, the range is found to be between 278 K and 
about 240 K. Above 278K, continuous Brownian diffusion dominates while below 240 K, hopping motion is found to 
control the mass transport jT]. Long time and long wavelength diffusion is in any event, insensitive to these details of 
the mass transport or relaxation mechanism found to be present in the crossover region but at intermediate to large 
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q, structural relaxation becomes sensitive to the gradually changing mechanism and thus provides a highly useful 
window to probe the details of the relaxation dynamics. Clearly, at very low temperature (much below MCT T° or 
onset temperature). and in the large q limit, relaxation should become independent of q if mass transport or relaxation 
is dominated by hopping. 

A characteristic feature of the dynamics in supercooled liquids is a pronounced non-exponential relaxation of the long 
time decay, commonly known as the a-relaxation. This relaxation is usually fitted to the Kohlrausch- William- Watts 
(KWW) function. 

<P( q ,t) = f( q )e-V^\ (1) 

where f(q),r(q), and (3(q) are the wavenumber dependent nonergodicity parameter, the relaxation time and the 
stretching parameter, respectively. This functional form exhibits highly non-trivial q dependence. When (f>(q, t) 
denotes the density correlator then f(q),r(q), and [3{q) oscillates in phase with the static structure factor, S(q). 
Although all the oscillations are in phase with S(q) each parameter exhibits different q dependence. We may also 
define another parameter T(q) which gives the combined effect of the q dependence of r(q) and /3(q), = T(q)~P( q \ 
The complete q dependence of the relaxation kernel is given by r(q). If T(q) has quadratic q dependence, </>(g, t) will 
be Gaussian. Weaker q dependence of T(q) leads to non-Gaussian <f>(q,t). 

In the simulation studies of Sciortino and coworkers on water ^ |S] and OTP [H], the variation of r(q) with q 
at different temperatures were studied. In their study in the large q limit, as the temperature is lowered, the q 
dependence of r(q) shows a crossover from quadratic to linear behaviour, 

r(q) cx q~ a . (2) 

At high T, a ~ 2 and its value decreases as the temperature is lowered. These authors also find that (3{q) varies 
from 1 at small q to 0.5 at large q. In their study, the q dependence of r(q) is always quadratic in the small q limit 
4 . while the dynamics in the large q limit and low temperature is stretched and non-diffusive 0J El E] • Similar 
subquadratic q dependence at low temperature has been observed in Brownian dynamic simulation of hard spheres 
[TJ. In a different simulation study on polymer melt the authors studied the q dependence of r(q) at three different 
temperatures [8]. In the large q limit, they tried to fit their data to mode coupling theory prediction, r(q) cx q~ x l h 
where b is the von-Schweilders exponent. They find that b — 0.75 which means that r(q) has a weaker than quadratic 
q dependence. However, if we analyze the results carefully we find that as the temperature is lowered there is further 
decrease in the exponent, i.e the q dependence of r(q) becomes even weaker. Thus all these studies clearly show a 
weakening of the q dependence of r{q) at low T. 
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According to the idealized mode coupling theory [T21 [TB] in the large q limit, 

r(q)-^ = T(q) oc q (3) 

IMCT also predicts that f3(q) varies from 1 at small q and then asymptotically saturates at about 0.5 at large q, thus 

Um^riq) ~ (r( 9 ))-V/*(«) « q- 1 /^) a g -a (4) 

TTms would mean that according to IMCT, the dynamics is diffusive even at large q. Although the timescale remains 
quadratic in q, due to the stretched dynamics, ((3(q) < 1) 4>(q,t) becomes non-Gaussian. 

The phenomenon of the weakening of the q dependence of r(q) as the temperature is lowered is often attributed 
to activated hopping dynamics |10j . However, the jump diffusion model which takes into consideration a distribution 
of jump length, predicts a much weaker q dependence in the large q limit and it also does not account for the non- 
trivial temperature dependence of the exponent a . So neither the continuous diffusive dynamics nor the activated 
dynamics by itself can explain the weakening of the q dependence as observed in simulations jH [5j El [7, 8 Both the 
theories fail to describe the progressive weakening of the q dependence of r(q) because in the supercooled liquid the 
dynamics is neither completely diffusive nor completely activated. As mentioned earlier, the dynamics changes rather 
smoothly from diffusive to activated and in the crossover region both mechanisms of relaxations are present pQ . 

In an earlier study |18j , we proposed a unified theoretical approach to combine ideal MCT with a nucleation based 
random first order transition theory treatment of the activated event. In a subsequent paper, [T] we showed that 
the unified theory can reproduce many of the observed experimental features, in addition to explaining the crossover 
behavior mentioned earlier. The theory predicts an interesting interplay between the two mechanisms, in the crossover 
regime. The earlier study [T] did not address the wave number dependence. In the present study we investigate the 
consequences of the theory for the wavenumber dependence of r(q) in the large q limit at different temperatures in 
the crossover regime. Here we will study the coherent intermediate scattering function. In the large q limit apart from 
the modulation of r(q) by the structure factor, S(q), the timescale of both the total (coherent) and the self (coherent) 
intermediate scattering function show similar q dependence. 

The combined theory suggests that an explanation of the wavenumber dependence of r(q) must essentially have 
contributions from both diffusive and activated dynamics as well as an interaction between these two different dy- 
namics which leads to hopping induced diffusion. The earlier study has shown that there is a gradual change in 
mechanism of relaxation from diffusive to activated, in the crossover regime. Here we show that change of mechanism 
of transportation or crossover has a wavenumber dependence. In the intermediate q regime the crossover is sharper 
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and takes place at a higher temperature. This q dependent crossover gives rise to the gradual weakening of the q 
dependence of the relaxation time. 



A. Theoretical scheme 

The previously developed theoretical scheme that combines the hopping (through RFOT theory) and continu- 
ous diffusive (described using MCT) dynamics can be applied to study the wavenumber dependence of the various 
dynamical quantities as described below. We also present a pure microscopic derivation of the memory function 

Let us say that an instanton pops up at a position R. The particles within the sphere of radius £ around the 
position R will be displaced by a Lindemann length, d/. Now the master equation describing the change in density 
p(r) due to a instanton popping up and due to the continuous diffusive dynamics together can be written as, 

6 P new (r, t + St) = Sp(r, t) + P hop (r) x Sty J dRd((r - R) < 

x [ dt' I dr' Sp(r,', t') x G(r, r', t - t') - 5p(r, t) 

J JV{R) 

-St { dt' { dr'8p(r,', t')K M c T (r, r', t - t') (5) 



where 6*((r — R) < £) provides the effect of the instanton felt at the position r, provided it is within the radius £. 
G(r,r't — t') is the Greens function which determines the effect of instanton in moving particles from position r to a 
new one r', typically a Lindemann length away, during time t and t + St. T>(R) determines the region where the effect 
of the instanton is felt [15] . Ph op (R) is the total rate of a hopping in the volume V. Kmct is the diffusive kernel 
which represents the diffusive dynamics as described by MCT [TJ [TS] . 
From Eq{5] the total scattering function can be written as [18] , 

<j)(q, z) = — — r (6) 

z + K R {q,z) 

where K]i(q,z) — K^^q^z) + Kmct(Qi z ) describes the renormalized kernel describing both the activated and the 
diffusive dynamics. However, in this form we cannot separately analyze the contribution from the diffusive and 
activated motion. Thus to describe the total intermediate scattering function we use an approximate form by making 
an approximation, 

^-cftMCTiqMhZtM ( 7 ) 

It has been demonstrated elsewhere that the equation of motion obtained for (f>(q, t) from cq|6] and eq|7] have very 
similar form and describe similar dynamics [18] . As mentioned earlier, the advantage of describing cf>(q,t) from eq{7] 
is that we can now separately analyze the activated and the MCT parts of the total dynamics. 



In describing the activated dynamics we consider that there is a distribution of the hopping barriers in the system 
arising from the entropy fluctuation |15j . Thus the total contribution from the multiple barrier hopping events is 
written as, 

^of C M) = / ^ hop ( q ,t)T staUc (AF)dAF 

= J e - tKh "^ AF) V staUc (AF)dAF 7 (8) 

where p statlc (^AF) is taken to be Gaussian, 

j,staUc( AF ^ _ 1 C -(AF-AF )/2(SAF) 2 /g\ 

yfMKAFf 

We call this distribution "the static barrier height distribution" . Here ^ op (i) describes the activated hopping dynamics 
for a single hopping barrier which has been derived in the earlier paper |18j, 

1 s + K hop (q, AFj 

where AF is the free energy barrier for hopping which determines the probability of a hop or the waiting time. 

K hop ( q ,AF) = ^^[l-G(q)] (11) 

Vp 

In the above expression of the hopping kernel, Phop is the average hopping rate which is a function of the free energy 
barrier height, AF and is given by P = -^exp(— AF/ksT) [H]. The free energy barrier is calculated from RFOT 
theory |16j . vq — §7r£ 3 is the region participating in hopping where the correlation length £ is calculated from RFOT 
theory. v p is the volume of a single particle in the system, is the Lindemann length. 

Supercooled liquid simulation studies have shown that there is in fact not a single length but rather a distribution 
of jump lengths which can be fitted to an exponential distribution |17j . The most probable jump length has been 
found to be of the order of Lindamann length with the value decreasing as the temperature is lowered. In this work 
we we consider the distribution of jump length, f (l) to be given by |llj . 

/„(/) = ll- 2 exp(-l/l ) (12) 

then the Greens function can be written as , 



1 

1 -fqHl 



l D is the most probable jump length which is found to be close to the Lindemann length. The hopping kernel can be 
written as 
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We now write the equation of motion for the MCT part of the intermediate scattering function, 4>MCT(q,t) which 
is now self consistently calculated with cj>(q,t), 

<t>MCT(q,t) + lq<t>MCT(q,t) + fi^0MCT (<?,*) 

+ / dt'n 2 q m 2 Jt')^MCT(q,t-t')=0 (15) 



where m q is the function of M variables cj>(q,t), (q=l,2,...M) and is given by, 

m,(t) = F q m)) 

= V(q;k,p)<p(k,t)<f>(p,t) 

k+p=q 

= [pS(q)S(k)S(p){q\kc(k)+pc(p)}} 2 /2q i }cj ) (k,t)<j ) (p,t) (16) 

k+p=q 
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T q is the mode coupling functional. The vertices V arc equilibrium quantities |19j . The double summation in Eq 
is approximated by a Riemann sum ["HJ] . where M = 100 and the cutoff wavevector q cuto ff a = 40.0. 

B. Results 

We first discuss the choice of parameters. Following our earlier work, we have taken the thermodynamic conditions 
appropriate for Salol, although the calculations are quite general. The system at high temperature is chosen to be just 
above the idealized MCT transition having T = 280K, p* = 0.99 and the lowest temperature T = 229K, p* = 1.27. 
The last value is suggested by the simulation studies [22] of deeply supercooled liquids with model Lennard- Jones 
potential, which shows that at p* = 1.27, T* = 0.5 the system is highly viscous. In van der Waals systems the large 
change in density drives the caging which can in principal arise not only due to change in density but due to transient 
bond formation which happens in many systems like water. Since we do not explicitly take into consideration the 
other sources of caging, so in a mean-field way, to mimic strong caging, the density is increased by a large amount. 
The density between T = 280JT and T = 22%K is calibrated in a linear fashion. With these values, we calculate the 
structure factor using the Percus-Yevik approximation with correction. The structure factor is used in the calculation 
of the microscopic MCT vertex. The parameters required to calculate the activated dynamics are discussed extensively 
in Reference [T] 

With the above mentioned parameters (j)(q,t) is calculated at different wavenumbers, q, and over the full time 
domain t. Numerical calculations were carried out using the scheme given in Ref. |20J. The longtime part of </>(g, t) at 



each wavenumber is fitted to the KWW function (eqjTJ, where all the KWW parameters, the form factor, f(q), the 
relaxation time, r(q) and the stretching parameter, (](q) are varied to get the fit. Similarly <pMCT{q^t) and 4>hop(lyt) 
are also fitted to KWW function to obtain the corresponding MCT and activated dynamics parameters. 

We plot logioT(q) against logio(q) at four different temperatures (figure 1). As the temperature is lowered, the 
wavenumber dependence of the relaxation time gradually becomes weak, similar to that observed in experiments 
[4j |5j |6j |7j [8]. In figure 2 we also plot log(T(q)q 2 ) vs. q to show the departure from an oscillatory quadratic 
wavenumber dependence to a sub-quadratic (almost linear) dependence at lower temperatures. The relaxation time 
obtained independently from the MCT and the activated parts are plotted in figure 3. In the large q regime, 
TMCT(q) as observed earlier in the case of idealized MCT [15] . shows a quadratic q dependence whereas Th op {q) is 
almost independent of q. The quadratic wavenumber dependence (r(q) oc 4y) is a signature of the continuous Brownian 
diffusion and the weak wavenumber dependence (r(q) oc is a signature of discontinuous activated hopping. Thus 
neither MCT nor activated dynamics by themselves would explain the subquadratic, almost linear, wavenumber 
dependence of the relaxation time. 

In an earlier study we have shown that in the crossover regime there is no abrupt change in the dynamics from 
diffusive to activated, the change is rather gradual p]. As we move closer to the laboratory glass transition tem- 
perature, the activated dynamics becomes more and more dominant. Here we investigate this crossover in the full 
wavenumber domain. The function At* = (tmct(?) — Ttotai(q))/T~totai(q) is plotted as a function of both temperature 
and wavenumber in figure 4. The higher the value of this function the stronger is the effect of hopping. We find that 
at high temperature, over the whole wavenumber regime the effect of activated dynamics is very weak. As the tem- 
perature is lowered, the effect of the activated dynamics shows strong wavenumber dependence. In the intermediate 
q regime the crossover is sharper and takes place at a comparatively higher temperature. In the large q regime the 
crossover is quite gradual. 

The wavenumber dependence when activated dynamics is included is largely a consequence of the wavenumber 
dependence of the hopping induced MCT relaxation time. An asymptotic analysis of the MCT dynamics provides 
the following expression for the timescale of hopping induced diffusive motion below the MCT transition temperature 

2T- 1 

-1 _ hop /-.-n 

t mct A/(g) 2 — 1 1 ' 

Here f(q) denotes the form factor or the Debye Waller factor (DWF). A denotes the coupling strength of the self 



coupling term and is a function of m q (t = 0) (given by eq. 16 1 . The above relation between hopping induced MCT 



relaxation time was derived at a single wavenumber q = 2n/a, without considering any coupling between wavenumbcrs. 
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The relation will be more complex for the present more complete treatment where all the wavenumbers are coupled. 
However it is easy to see that even in the decoupled limit, the relaxation time becomes longer as the DWF increases. 
Thus at large q although Thop is independent of q, the hopping induced diffusive dynamics has a q dependence and 
the relaxation time is longer for wavenumbers where the DWF is larger. 

In figure 5 the total relaxation time, r tota /, the hopping induced MCT relaxation time, tmct an d the relaxation 
time of the activated hopping dynamics, Thop are plotted at three different temperatures. At high T, for all q, 
as discussed earlier, the dynamics is primarily determined by MCT. However, at low T and in the intermediate 
wavenumber regime the MCT relaxation time becomes very slow (due to caging effect) thus the relaxation primarily 
occurs via activated dynamics. As observed in figure 5, in the large wavenumber regime, even at low T both MCT 
and activated dynamics have competing timescales and the relaxation has contributions from both. This wavenumber 
dependence of the effect of hopping and stronger effect of activated dynamics at intermediate wavenumbers would 
imply that the growth of relaxation time with lowering of q will not be as strong (quadratic) as predicted by the MCT 
dynamics. 

This theory also predicts that even in the supercooled regime there is a limited range of q where non-Gaussian 
behaviour can be observed. This range, which is primarily positioned in the intermediate wavenumber regime, 
becomes asymmetrically wider as the temperature is lowered. As observed in figures 4 and 5, the effects of hopping 
dominated regime extends to higher q values with lowering of T. Thus the non-Gaussian regime also extends to higher 
values of q. However in the lower wavenumber side even when the hopping dominated regime is reached the non- 
Gaussian wavenumber dependence will be restricted only above q * a ~ 2n. This is because, the activated dynamics 
itself predicts a Gaussian wavenumber dependence at small q. 

C. Discussion 

The unusual wavenumber dependence of the relaxation time observed in supercooled liquids has heretofore remained 
largely unexplained |U El EJ [3 |S] . This article addresses this problem by extending the bridged theory [TJ [TH] to the 
full wavevector plane. We show that the full wavevector dependent unified theory explains the gradual weakening of 
wavenumber dependence of the relaxation time as the temperature is lowered. The theory predicts that for systems 
where the caging effect grows strongly with lowering of temperature, the MCT relaxation time in the intermediate 
wavevector regime slows down faster and the crossover to the activated dynamics dominated regime happens earlier. 
We have shown that this crossover to the purely activated regime is strongly wavenumber and temperature dependent. 
The crossover to the activated regime happens earlier at intermediate wavenumbers, whereas at the same temperature 
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and at larger wavenumbers, the MCT contribution still remains substantial. Because of this wavenumber dependent 
crossover the relaxation time at low T is governed mostly by activated dynamics at intermediate q and both by 
activated and MCT dynamics at larger q. Hence the relaxation time does not grow as strongly with lowering of q as 
expected from only the diffusive MCT dynamics. The theory also predicts that if the caging effect is not so strong 
so as to significantly slow down the MCT dynamics (compared to hopping) we should observe quadratic wavenumber 
dependence even at lower temperature. The analysis also reveals that for systems where the caging effect is stronger, 
the transition from MCT dominated regime to the activated dynamics dominated regime should be sharper. This 
would mean that the phenomenological MCT transition temperature, T/ lt will appear to be closer to the microscopic 
MCT transition temperature T®. In other cases where the caging effect is not very strong or does not grow very fast 
upon cooling the transition takes place over a wider temperature regime and the separation between T[ %t and T° 
should be wider. 

It is now believed that supercooled liquids are characterized by a growing dynamic correlation length, £. For 
activated transitions in pure RFOT theory, the growth of this correlation length is given by the configuration entropy, 
£ cx s c 2 ^ 3 16J. At high temperature where this dynamic correlation length is small, relaxation is described by the 
MCT. As the temperature is lowered, £ grows, and as discussed here, with the lowering of T, activated dynamics 
also becomes progressively a major contributor to the total relaxation dynamics. It is thus reasonable to ask if 
there is a correlation between £ and the non-Gaussian wavenumber dependence. The present analysis suggests that 
while a growing correlation length is essential to describe the timescale of activated dynamics and thus the non- 
Gaussian dependence is related to £, there appears to be no one-to-one relation between £ and the non-Gaussian 
dependence for the following reason. The non-Gaussian dependence is not described by activated dynamics alone but 
it is MCT and activated dynamics and their interplay which gives rise to the non-Gaussian wavenumber dependence 
of the relaxation time. Since MCT dynamics is independent of £ and depends on many other factors like caging, the 
relationship between the correlation length £ and the non-Gaussian wavenumber dependence may not be universal. 
For the system we have studied, it is fair to say that sub-quadratic wavenumber dependence becomes significant only 
when the correlation length £ > 2. 
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FIG. 1: The a relaxation timescale r(g) plotted as a function of q at different densities and temperatures. The r(q) values are 
scaled such that at q = 8.6 they have similar values. r(q) shows a weaker q dependence as the temperature is lowered. 
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FIG. 2: r(q)q 2 plotted as a function of q at different densities and temperatures. The r(q)q 2 values are scaled such that at 
q — 12.5 they have similar values. As the temperature is lowered, r(q) shows a crossover from quadratic to sub-quadratic 
(almost linear) q dependence. 




FIG. 3: The a relaxation timescale r(q) for the MCT part, Tmct an d the activated part, Th op are plotted as a function of q. 
tmct shows a strong variation with q whereas Thop has a much weaker q dependence. The relaxation times are calculated at 
T = 2QZK and p = 1.16 
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FIG. 4: Wavenumber dependent crossover to the activated regime. The function At* = (tmct — T tota i) I Tmct is plotted 
against both temperature and wavenumber. The closer the value of the function to unity the more is the effect of activated 
dynamics. In the intermediate wavenumbers, the crossover to the activated dynamics dominated regime is sharper and takes 
place at a higher temperature. In the large q regime this crossover takes place over a wide temperature range. 
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FIG. 5: (a)The total relaxation time, T tota i, the hopping induced MCT relaxation time, tmct and the relaxation time of the 
activated hopping dynamics, r; lop are plotted against time at T=280K. The total relaxation is primarily determined by the 
hopping induced diffusive dynamics. (b)Same as in (a) but at T=263K. The total relaxation in the intermediate wavenumber is 
primarily determined by the hopping dynamics whereas in the long q regime the total dynamics follows the relaxation timescale 
of the hopping induced diffusive dynamics. (c)Same as in (a) but at T=247K. Here too the total relaxation in the intermediate 
wavenumber is primarily determined by the hopping dynamics and only in the very large wavevector regime MCT plays a 
dominant role. 



